On the ground state energy scaling in quasi-rung-dimerized spin ladders 



(N 

o 

< 
O 



u 

-f— > 
-4-* 



I 

c 

o 
o 



>: 

(N, 
OV 

06 

o: 



P. N. Bibikov 1 

1 Saint-Petersburg State University 
(Dated: August 31, 2012) 

On the basis of periodic boundary conditions we study perturbatively a large N asymptotics (N 
is the number of rungs) for the ground state energy density and gas parameter of a spin ladder with 
slightly destroyed rung-dimerization. Exactly rung-dimerized spin ladder is treated as the reference 
model. Explicit perturbative formulas are obtained for three special classes of spin ladders. 
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I. INTRODUCTION 

Phase structure of frustrated spin ladders and spin lad- 
ders with four-spin terms has been intensively studied 
in the last decade both theoretically and numerically-ir— . 
Among other phases the mathematically most simple one 
and at the same time, probably, the one most interest- 
ing for physical applications is the so called rung-singlet 
(or rung-dimerized) phased. Within it the ground state 
may be well approximated by an infinite tensor product 
of rung-dimers (singlet pairs) 



\0)r-d = ®n|0), 



(1) 



This state will be an exact ground state only for rather 
big antiferromagnetic rung-coupling and under a special 
condition on the coupling constants^. The latter has no 
physical background and thus there are absolutely no 
grounds to assume its relevance for real compounds. Nev- 
ertheless it is a common opinion that for rather big anti- 
ferromagnetic rung coupling a spin ladder should still re- 
main in the rung-singlet phase. This means that all phys- 
ical properties of such a ladder may be obtained pertur- 
batively on the basis of the "bare" ground state (1) and 
its excitations. Together with verification by machinery 
calculations this approach should give a comprehensive 
description of the rung-singlet phase. A machinery calcu- 
lation will provide excellent tests for suggested formulas 
while a perturbative formula will give a right direction 
for numerical research and interpretation of the obtained 
data. 

Such approach has two main difficulties. First of all 
a general spin ladder model is non-integrable and al- 
though one- and two-magnon states may be readily de- 
rived within Bethe Ansatz, three-magnon states are ob- 
tained now only for five special integrable models^. The 
second difficulty originates from the fact that an analyti- 
cal result is usually obtained for infinite ladder however in 
a numerical calculation a ladder has a finite size. Hence 
in order to use a mashinery calculation for verification of 
an analytical result one havs to perform a correct extrap- 
olation of the numerical data. This means that utilizing a 
finite number of numerical estimations Jn of some value 
/ (N the number rungs of the ladder) it is nesessary to 
estimate the limit foo — limAr^oo /jy. On this way, in 
addition to a number of sequence transformation meth- 



ods improving the convergence^, one has to be guided by 
some extrapolation formula. The latter may be guessed 
by an analysis of numerical data^, or suggested theoreti- 
callly on the basis of conformal field theory 10 predictions, 
or on some other argumentation^. 

Taking an exact rung-dimerized spin ladder as a ref- 
erence model, it is natural to treat the ground state of 
a spin ladder with violated rung-dimerization as a dilute 
magnon gasi^. Its consentration (gas parameter) 



= lim 

JV-s-oc 



Pn, Pn = 



(0JQJ0) 
N 



(2) 



(Q is a magnon number operator (13)) and energy density 



E = Eoo = lim E N , 

N-s-oo 



En — 



|g|0) 
N ' 



(3) 



turns to zero for an exact rung-dimerized spin ladder 
and hence they should be good governing parameters for 
a perturbation theory based on the gas approximation. 
Perturbative expressions for p and E were derived in Ref. 
12. In the present paper assuming periodic boundary con- 
ditions we obtain in three special cases the corresponding 
extrapolation formulas for pn and En- 
The two formulas 



E 



N 



E x + (-1) N A 



,-N/N 



E 



N 



A_ 



(4) 
(5) 



(A and Nq are free parameters) have already been 
suggested correspondingly for ope n 13 ' 14 and periodic)^ 
boundary conditions. The expression (4) was implied ad 
hoc, while Eq. (5) follows from conformal theory argu- 
mentation. The perturbative formulas obtained below 
for three special classes of spin ladders have a rather dif- 
ferent form 



E N = E 00 + {A+{~l) N B)e 



-{N-X)/N 



(6) 



II. DESCRIPTION OF THE MODEL 

We shall use an equivalent representation^^- 
H = H + J 6 V, 



(7) 
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of the spin ladder Hamiltoniani — . Here Jq is a pertur- 
bation parameter and 



N 



H = ^ J lQn + J 2 (*n •*„+!+*„• 

n=l 

+ J^QnQn+1 + J 4 + J 5 (s 

JV 

V = ^K, n+1 , (8) 



1 



Sri — Sl )TJ + S 2jJl , Q n — 2^"' 



(9) 



(Si tn for i = 1,2 are spin-1/2 operators associated with 
ri-th rung). The local operators 

= 2 (^1,71 ~ ^2,n) — i[Sl,n X S 2j „], 

*„ = ^(Si,„ - S 2 ,„) + i[Si,„ x S 2 , n ], (10) 

may be interpreted as (neither Bose nor Fermi) creation- 
annihilation operators for rung-triplets. Namely 

#-io)„ = ii>-, n\it = o, 

*£|0>„ = 0, #°\l) b n = S ab \0) n . (11) 

From (8) and (9) readily follows^ that 

[H ,Q]=0, (12) 

where the operator 

Q=5ZOn, (13) 



according to relations 



(14) 



has a sence of the number operator for rung-triplets^. 

For rather big J\ (for example nesessary should be^& 
J\ > J2) vector (1) is the zero energy (i?o|0) r _d = 0) 
ground state for Hq, whose physical Hilbert space splits 
into a direct surn^&i^ 



(15) 



The subspace H° is generated by the single vector (1). 
According to (2), (3) and (8) 



Pn 



dE 



N 



(16) 



Since V : |0) r _c; — > H 2 , a perturbative treatment of 
the term J 6 V gives 



N 



\{»\V\0)r-d\ 

E(fj,) 



+ o(4), (17) 



where all the states in the sum have zero total spin 
and quasimomentum. In the N — > 00 limits 

g 00 = -e(Ag-i) 3 5 ( ^"~ 1) 



3J| 



1 



J 2 |A^-l| + 2Ao V /J 1 2 - Jf 
[2A Ji + (Aq + 1) J 2 ] 



(18) 



where 0(x) = 1 for a; > and 0(x) = for x < and 
J 3 - 2 J 4 + 4 J 5 



A ° = 2~T 2 

Ebound = 4Ji + 2 J 2 



(19) 
(20) 



III. A FINITE- iV TWO-PARTICLE PROBLEM 

A zero total spin and quasimomentum two-magnon 
state has the following general form, 

|2 -magn)= ]T a{n - m)...\l) a m ...\l) a n ... (21) 

l<m<n<N 

The dimension of the corresponding Hilbert space is N/2 
for even N and (N — l)/2 for odd. The wave function 
a(n) should be normalised 

JV-l 

£ (iV - n)\a(n)\ 2 = ]T |o(n - m)\ 2 = -, (22) 

n— 1 m<n 

and satisfy the periodicity condition a(n — m) = a(m + 
N — n) or shortly 



a(n) = a(N — n). 



(23) 



Performing a substitution n — > jV — rt and using (23) one 
can obtain from (22) 



JV-l 



JV— 1 



J^n\a(N-n)\ 2 = £ n|a(n)| a = ± (24) 



n— 1 n— 1 

Together (22) and (24) result in 

JV-l 



n=l 



3iV' 



(25) 



The Schrodinger equation gives 
4Ji<z(n) + 2J 2 [a(n - 1) + a(n + 1)] = Ea(n) (26) 
for 1 < n < N — 1 and 

2(2 Ji + J 2 A )o(l) + 2J 2 a(2) = £Ja(l), (27) 
for n = 1. 
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General solution of the system (26), (27) has the form 
1 



a(n, z) = 



sfW) 
and dispersion 

E{z)=U l +2J 2 (z + -^j. 



(28) 



(29) 



The normalization constant Z(z) ensures condition (25). 
The parameter z corresponds to relative quasimomentum 
of magnon pair and satisfy an equation 



JV-l 



Apz-l 
z- A 



Ap-l/z 
A -z 



(30) 



The latter is invariant under complex conjugation and a 
duality symmetry 



(31) 



which according to (28) is related to multiplication of 
the wave function on (-1). Hence for even N the roots of 
(30) are joined in dual pairs, while for odd N there is an 
additional autodual root z = — 1. 

In the three special cases A = — 1, A = 1 and A = 
Eq. (30) may be solved explicitly. Denoting the corre- 
sponding solutions RS 11 j , Vj and Wj respectively one has 



_ „(2j+l)i7T/(JV-l) 



Uj = c 



S 2jiw/(N-1) 

J2j+1)17T/N 



0,...,N- 
0,...,N- 2, 
0,...,iV-l, 



2, (A = -l), 
(A = 1), 
(A = 0). (32) 



Taking into account that all the roots (32) lie in a unite 
circle one may readily get 



Z(z) = 3JV(JV-l)(l-A 2)(l--^), A = ±l, 
Z(z) = 3N 2 , A = 0, (33) 

and then 



\a(n,z)\ 2 = 



3N(N - 1) 
A = ±1, 



2 + Ao(, 2 - 1 + 



M">*)\' = 3^( 2 -- 2 "-^r)' A " = 0. (34) 



IV. EXACT RESULTS AT A = AND A = ±1 

Let \z) be the state related to wave function (28). From 
(9) and (21) follows that 



(z|f|0) r ,_ d | 2 =9iV 2 |a(l,z)| 2 . 



(35) 



For the evaluation of En one has to perform in (17) a 
summation over all duality pairs of roots. Since both the 
roots in a pair give the same contribution this is equiv- 
alent to inserting the factor 1/2 before summation over 
all roots. Hence (17) and (35) result in 



£jv(A ) 



--JlG N (A ) + o(Ji), 



(36) 



where 



Gjv(-I) 



G N (1) 



G N (0) 



1 ^ 2 - ( Uj ; + 1/Uj) 



JV-2 



N 



-Y — 

1 ±* 2Jj 



+ J 2 (u j + l/u j ) J 2 (N-1) * 



+ J 2 ( J- ' J+ ' 

i ^ 2 + ( Vj + i/ Vj ) l ^ r 

N — 1 4-t 2Ji + Joivi + 1/v-i) MN-D^A 



N-l 



J 



N-l 



JV-2 



Jl + J 2 



- J2 I 

~~J 2 \ 



y/J't - J 2 V J- - u 3 J + - u 3 



o 2Ji + J 2 (ty + Ja(JV-l) . „ 



1 - 



J1-J2 



l_\ J1-J2 



■h 



yJ\ - J 2 



J 



N-l 



J ±l ) 

l^N—l jN-1 + (_^N-l) 



T N-1 



N-l 



N E Wj [J 2 w] - 2J x w 3 + J 2 ) J\N 



3=0 



Ji \J J\ — J 2 ( J 



E f* - #h + -) - \1^~JI(t ± ~r~^~ — ) 

j-^l 2 V J Wj J V 1 2 \ j_ - Wj J + - Wj ) 



J 2 

J 2 



T 2 

J 2 



, JN jN 

b^+T ~Jf + lJ 



(37) 



and 



J± = 



-J 1 ± y /j^J i 



(38) 



In (37) we used for calculations the formulas 

N-2 , f TtT -.wjV-2 



E 

3=0 
N-2 



1 



J 



E T _ „. 



(TV - 1) J 

J^- 1 + 1 ' 

(AT - 1) J^" 2 
fiv-i Z uw ; 



4 



which may be proved according to the following argu- 
mentation. The sums in (39) are fractions whose numer- 
ator and denominator are symmetric polynomials with 
respect to Uj, vj and Wj respectively. However according 
to (30) all these polynomials exept 



U --U N -2 
W ...W N -! = 

are equal to zero. 
From equality J + J_ 



(-1)"- 



Vq—Vn-2 



= 1, 



(40) 



1 readily follows 



J 



JV-l 



'+1 

jN-1 



jN-1 
J?- 1 + 1 



1 - J 



JV-l 



1 + J^- 1 ' 



r JV-l 
J + 



JN-1 + (_1)JV-1 + (-l)JV-l 

1 - (- J+)"- 1 



(41) 



l + (-J + )iv-i' 
Using (41) one may readily reduce Eqs. (37) to the form 



Gjv(-I) 



1 

^2 



■A — ^2 



G N (1) = - 



1 + 

1 



JV-l 



1 



(-./+) 



JV-l 



Gjv(0) 



2 

J2 



Jl 

J 2 



Ji-J: 

Ji + J 2 ' 1 + (- J+) N - 



i 



N 



(42) 



It may be readily observed that the corresponding val- 
ues for E oo (A ) agree with Eq. (18). The scaling law 
has the form (6) with 



= 0, B(-l) 



3J| / Ji + J<i 



3J fi J\ — J2 



2J 2 \ J1 + J2 



2J2 V J\ _ J2 ' 

5(1) =0, 



3J, 



A(0) = 0, B{0) = -§-Jj*-J*, 

J 2 



at J2 > and 



^ ^ 2J2 \ Ji — J2 



3 J| / Ji + J2 



, B(-1)=0, 



= 0, 5(1) = -^ 



3 Jr Ji — J2 



2J 2 \ Ji + J 2 ' 



N 



ln|J 2 | -ln(Jx - yO? - J|) 



(43) 



3 7 2 / 

^(0) = ^V J i- J 2' B (°) = °> ( 44 ) 



at J2 < 0. In both the cases 



(45) 



The corresponding formulas for p/y have the similar 
form and may be readily obtained from (16). 
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